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ABSTRACT 

The transfer of polarized radiation in stochastic synchrotron sources is explored by means of analytic 
treatment and Monte Carlo simulations. We argue that the main mechanism responsible for the circular 
polarization properties of compact synchrotron sources is likely to be Faraday conversion and that, 
contrary to common expectation, a signihcant rate of Faraday rotation does not necessarily imply strong 
depolarization. The long-term persistence of the sign of circular polarization, observed in many sources, 
is most likely due to a small net magnetic flux generated in the central engine, carried along the jet axis 
and superimposed on a highly turbulent magnetic held. We show that the mean levels of circular and 
linear polarizations depend on the number of held reversals along the line of sight and that the gradient 
in Faraday rotation across turbulent regions can lead to “correlation depolarization”. Our model is 
potentially applicable to a wide range of synchrotron sources. In particular, we demonstrate how our 
model can naturally explain the excess of circular over linear polarization in the Galactic Center and the 
nearby spiral galaxy M81 and discuss its application to the quasar 3C 279, the intraday variable blazar 
PKS 1519-273 and the X-ray binary SS 433. 

Subject headings: polarization: circular - stars: individual (SS433) - Galaxy: center - galaxies: 

individual (M81) - quasars: individual (3C 279) - BL Lacertae objects: individual 
PKS 1519-273 


1. INTRODUCTION 

Polarization has proven to be an important tool in AGN 
research. In principle, linear and particularly circular po¬ 
larization observations of synchrotron radiation may per¬ 
mit measurements of various properties of jets such as: 
magnetic held strength and topology, the net magnetic 
hux carried by jets (and hence generated in the central 
engine), the energy spectrum of radiating particles, and 
the jet composition (i.e., whether jets are mainly com¬ 
posed of e+ — e~ pairs or electron-proton plasma). The 
renewed interest in polarization of compact radio sources 
stems from two recent developments. First, Bower et al. 
(1999) detected circular polarization using the Very Large 
Array (VLA) in the best supermassive black hole candi¬ 
date, the Galactic center (Sgr A*). This discovery was 
quickly conhrmed by Sault and Macquart (1999) using 
the Australia Telescope Compact Array (ATCA). Circu¬ 
lar polarization was also detected in the celebrated X-ray 
binary system SS 433 (Fender et al. 2000). Moreover, 
the Very Long Baseline Array (VLBA) has now detected 
circular polarization in as many as 20 AGN (Wardle et al. 
1998; Homan and Wardle 1999). Second, it is now possible 
to measure circular polarization with unprecedented accu¬ 
racy of 0.01% using the ATCA (Rayner et al. 2000). This 
dramatic improvement in the observational status of polar¬ 
ization measurements has also brought new questions. For 
example, there is now growing observational evidence that 
the sign of circular polarization is persistent over decades 
(Komesaroff et al. 1984; Homan and Wardle 1999), which 
indicates that it is a fundamental property of jets. An¬ 
other problem, which has not been satisfactorily explained 
as yet, is how to reconcile the high level of circular polar¬ 


ization with the lower value of linear polarization in Sgr A* 
(Bower et al. 1999) and M81* (Brunthaler et al. 2001). 
Indeed, there is not even a general consensus on the mech¬ 
anism responsible for the circular polarization properties 
of jets (Wardle et al. 1998). 

In this paper we attempt to solve some of the theoret¬ 
ical puzzles. The paper is organized as follows. In the 
next section we summarize the most important observa¬ 
tional facts. In Section 3 we briefly discuss mechanisms 
for producing circular polarization and argue that the most 
likely process is “Faraday conversion”. Section 4 presents 
our model for polarization and in the subsequent sections 
we compare Monte Carlo simulations with analytic results 
and discuss general implications for observations as well 
as specific observational cases. We summarize our conclu¬ 
sions in Section 7. 

2. OBSERVATIONAL TRENDS 

Compact radio sources typically show a linear polariza¬ 
tion (LP) of a few percent of the total intensity (Jones et al. 
1985). This is much less than the theoretical maximum for 
synchrotron sources, which can approach 70% in homoge¬ 
neous sources with unidirectional magnetic field. There¬ 
fore, magnetic fields in radio sources are believed to be 
highly inhomogeneous, although the nonvanishing linear 
polarization is in itself an indirect indication of a certain 
degree of ordering of the held. Although the precise topol¬ 
ogy of the magnetic held in jets is not known, there are 
other compelling reasons, both theoretical and observa¬ 
tional, to believe that magnetic helds are indeed partially 
ordered. On the observational side, these conclusions are 
based on measurements of the orientation of linear polar- 
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ization (E vector), which reveal coherent structures across 
jet images. This indicates that magnetic fields, which are 
predominantly perpendicular to the electric vectors, are 
also preferentially aligned, although in different sources 
or in different parts of the jet the magnetic fields can be 
mainly orthogonal (Wardle et al. 1998) or parallel (Jones 
et al. 1985; Rusk and Seaquist 1985) to the projected 
jet orientation. From the theoretical point of view, ordered 
jet magnetic field is expected when shocks compress an ini¬ 
tially random field (Laing 1980, 1981; Marscher and Gear 
1985; Hughes et al. 1989; Wardle and Roberts 1994) (B 
perpendicular to the jet axis) or when such initial fields 
are sheared to lie in a plane (Laing 1980, 1981; Begelman, 
Blandford and Rees 1984) (B parallel to the jet axis). 

Circular polarization (CP) is a common feature of 
quasars and blazars (Rayner et al. 2000; Homan et al. 
2001 ), is usually characterized by an approximately flat 
spectrum, and is generated near synchrotron self-absorbed 
jet cores (Homan and Wardle 1999). CP is detected in 
about 30%-50% of these objects. Measured degrees of CP 
are generally lower than the levels of linear polarization 
and usually range between 0.1 and 0.5% (Homan and War¬ 
dle 1999; Homan et al. 2001). As reported by Macquart 
et al. (2000) for the intraday variable source PKS 1519- 
273, the CP of variable components of intensity can be 
much higher than the overall circular polarization levels. 
Observations of proper motion of CP-producing regions 
in the quasar 3C 273 (Homan and Wardle 1999) sug¬ 
gest that circular polarization is intrinsic to the source, as 
opposed to being due to foreground effects. Most impor¬ 
tantly, comparisons of CP measurements made within the 
last 30 years (Weiler and de Pater 1983; Komesaroff et 
al. 1984; Homan and Wardle 1999) with the most recent 
observations reveal that, despite CP variability, its sign is 
a persistent feature of AGN, which must therefore be re¬ 
lated to a small net unidirectional component of magnetic 
field in jets. 

3. MECHANISMS FOR PRODUCING CIRCULAR 
POLARIZATION 

The most obvious candidate for explaining circular po¬ 
larization of compact radio sources is intrinsic emission 
(Legg and Westfold 1968). Although intrinsic CP is 
roughly TTc^mt ~ where 7 is the Lorentz factor of ra¬ 
diating electrons, in a realistic source it will most likely 
be strongly suppressed by the tangled magnetic field and 
possibly the emissivity from — e~ pairs, which do not 
contribute CP. Specifically, 7rc,int ~ 7 “H'S«/^rms)/pair < 
1 %, where and are the unidirectional component 
of the magnetic field projected onto the line-of-sight and 
the fluctuating component of the field, respectively, and 
/pair = {n~ —n~^)/{n~ +n~^) < 1. Therefore, intrinsic CP 
appears to be inadequate to explain the observed polariza¬ 
tion. Other mechanisms have also been proposed, among 
which the most popular ones are coherent radiation pro¬ 
cesses (Benford and Tzach 2000), scintillation (Macquart 
and Melrose 2000) and Faraday conversion (Pacholczyk 
and Swihart 1975; Jones and O’Dell 1977a; Jones 1988; 
Wardle et al. 1998). The first of these mechanisms pro¬ 
duces polarization in a narrow frequency range which now 
seems to be ruled out by multiband observations. The re¬ 
cently proposed scintillation mechanism, in which circular 


polarization is stochastically produced by a birefringent 
screen located between the jet and the observer, fails to 
explain the persistent sign of circular polarization as the 
time-averaged CP signal is predicted to vanish. The last 
mechanism — Faraday conversion — seems to be the most 
promising one and in the next subsection we discuss it in 
more detail. 

3.1. Faraday rotation and conversion 

The polarization of radiation changes as it propagates 
through any medium in which modes are characterized by 
different plasma speeds. In the case of cold plasma the 
modes are circularly polarized. The left and right circu¬ 
lar modes have different phase velocities and therefore the 
linear polarization vector of the propagating radiation ro¬ 
tates. This effect is called Faraday rotation and it is often 
used to estimate magnetic field strength in the interstellar 
medium and to estimate pulsar distances. Note that Fara¬ 
day rotation does not alter the degree of circular polar¬ 
ization, since any circular polarization can be decomposed 
into two independent linearly polarized waves. Faraday ro¬ 
tation is a specific example of a more general phenomenon 
called birefringence. In a medium whose natural modes are 
linearly or elliptically polarized, such as a plasma of rel¬ 
ativistic particles, birefringence leads to the partial cyclic 
conversion between linearly and circularly polarized radi¬ 
ation as the phase relationships between the modes along 
the ray change with position. This effect is best visualized 
by means of the Poincare sphere (Melrose and McPredhan 
1991; Kennet and Melrose 1998). An arbitrary elliptical 
polarization can be represented by a vector P with its tip 
lying on the Poincare sphere and characterized by Carte¬ 
sian coordinates (Q, U, V)/I, where Q, U, V and / are the 
Stokes parameters (see Fig. 1). Thus, the north and south 
poles correspond to right and left circular polarizations 
and points on the equator to linear polarization. Different 
azimuthal positions on the sphere correspond to different 
orientations of the polarization ellipses. The polarization 
of natural modes of the medium is represented by a di¬ 
agonal axis, whose polar angle measured from the vertical 
axis depends on whether the medium is dominated by cold 
(0°) or highly relativistic particles (90°). As the radiation 
passes through the medium, birefringence causes the tip 
of the polarization vector to rotate at a constant latitude 
around the axis of the natural plasma modes. In this pic¬ 
ture, Faraday rotation corresponds to the case where the 
natural modes axis is vertical and the polarization vector 
P rotates around it. Note that, even if radiation initially 
has no circular polarization (i.e., P lies in the equatorial 
plane) and then encounters a medium in which the normal 
modes are elliptical, it will develop an elliptically polar¬ 
ized component. An interesting property of a relativistic 
birefringent plasma is that it can generate circular polar¬ 
ization even if it is composed almost entirely of electron- 
positron pairs. At first this may seem paradoxical, as one 
would expect electrons and positrons to contribute to CP 
with opposite signs. However, despite the fact that the 
intrinsic CP in such a case is indeed close to zero, some 
additional CP can be produced by conversion of linear po¬ 
larization (Sazonov 1969; Noerdlinger 1978). In terms 
of the Poincare sphere, this situation corresponds to the 
normal modes axis pointing close to but not exactly in the 
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equatorial plane. Therefore some conversion of intrinsic Averaging the transfer equations over orientations of 

linear polarization may occur provided that there is some the magnetic field, we obtain the following asymptotic 
imbalance in the number of electrons and positrons. expressions for large synchrotron optical depth: 


3.1.1. Strong rotativity limit 

Strong departures from mode circularity occur only 
when radiation propagates within a small angle ~ vl/i' 
of the direction perpendicular to the magnetic field, where 
vl = eB 12'KmeC. Therefore radiative transfer is often per¬ 
formed in the quasi-longitudinal (QL) approximation. If 
the normal modes are highly elliptical then the opposite, 
quasi-transverse (QT), limit applies (Ginzburg 1961). In a 
typical observational situation it is usually assumed that 
Faraday rotation within the source cannot be too large, 
as this will lead to the suppression of linear polarization. 
However, this constraint does not prevent rotativity from 
achieving large values locally as long as the mean rotativ¬ 
ity, i.e., averaged over all directions of magnetic field along 
the line of sight, is indeed relatively small. Such a situa¬ 
tion may happen in a turbulent plasma. Some effects of 
turbulence on polarization were discussed by Jones (1988) 
who neglected a uniform magnetic field component and by 
Wardle et al. (1998), who presented results for the case 
of a small synchrotron depth. Technically, the strong ro¬ 
tativity regime is equivalent to the QL limit and in this 
paper we build our model on this approximation. 

4. MODEL FOR POLARIZATION 


T+(2p-l)C,Q = J (1) 

Q + {CU) + {2p-l)(:<ri = {2p-l)egJ ( 2 ) 

I7-(C:Q) + (2p-1)C;F = 0 (3) 

F-(2p-l)C;c7=0. (4) 

Note that <C* is not statistically independent from U and 
Q due to the gradient in Faraday rotation across each cell. 
The correlation in eq. (2) then reads: 

{CU) = 5(U + {CU), (5) 


where Q and U denotes the fluctuating part of Q and 
U. An analogous relation holds for Q in eq. (3). We ne¬ 
glect the term (Q U) in eq. (4) as convertibility is a much 
weaker function of plasma parameters than rotativity. It 
will be shown in Section 4.1.3 that the correlations {QU) 
and (QQ) tend to zero as the number of field reversals 
along the line of sight increases and that 

-{CQ)/U = {CU)/Q = ^. (6) 


We consider a highly tangled magnetic field with a very 
small mean component which is required to determine the 
sign of circular polarization. From a theoretical view¬ 
point, we would expect some net poloidal magnetic field, 
either originating from the central black hole or from the 
accretion disk, to be aligned preferentially along the jet 
axis. Specifically, from equipartition and flux freezing ar¬ 
guments applied to a conical jet (Blandford and Konigl 
1979) we get ~ ~ firms oc r~^ where 

r is the distance along the synchrotron emitting source 
and the symbols || and T refer to magnetic fields paral¬ 
lel and perpendicular to the jet axis, respectively. From 
the flux-freezing argument applied to the small paral¬ 
lel bias in the magnetic field we obtain {B±) ~ 0 and 
(fi||) oc oc Jfirms, where 5 = fiu/firms <C 1 is the ratio 
of the uniform and fluctuating components of the magnetic 
field. 

4.1. Mean Stokes parameters in the presence of field 
reversals 

We solve the radiative transfer of polarized radiation in 
a turbulent plasma by adopting transfer equations for a 
piecewise homogeneous medium with a weakly anisotropic 
dielectric tensor (Sazonov 1969; Jones and O’Dell 1977a). 
Details of the transfer equations are given in the Ap¬ 
pendix. We assume that the mean rotativity per unit syn¬ 
chrotron optical depth (Q) = Sf and that (sin2^) = 0 
and {cos2(j)) = 2p — 1, where 0 < p < 1 is a parameter 
describing the polarization direction and degree of order 
in the field. We also assume that circular absorptivity (v 
and circular emissivity e„ are both negligible. 

4.1.1. Large synchrotron depth limit 


This implies that the mean levels of circular and lin¬ 
ear polarizations will also depend on the number of the 
field reversals along the line of sight. Setting Wq = Q/I, 
¥„ = {/// and Tty = V //, we get mean normalized Stokes 
parameters from equations (l)-(4): 

_ *9 

_ _ (2p -l)[l-^ + {2p- ifiCg ](e, - Cq) 

_ _ (2p- l)JC(e9 - Cg) 

V 

_ {2p-lfXtK{eq-Cq) 

where 

= 1 + 2^ + (2p - l)2(Cf - C,e,) + J'C" + 

[^-{2p-lfC,qeq][i+{2p-l)%^], (10) 

Note that 

^ = (2p-i)c; (11) 

T^u 


(7) 

( 8 ) 
(9) 


^ ^ (2p - 1)C/C ^^2) 

^q l + f + {2p- lYl*q ' 

Therefore, circular polarization Tty = —Ity can domi¬ 
nate over linear polarization tT; = (if^ il '^C ^ 

{2p — 1)Q > 1, where the last inequality holds for a large 
number of field reversals along the line of sight. 
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4.1.2. Small synchrotron depth limit 

In the limit of small synchrotron depth and for > 1, 
we approximately have: 


dl 

d^~'^ 

^+SCU+{CU) = {2p-l)egJ 

^ - iCQ) = 0 

^-(2p-l)C;[7=0. 


Introducing X = Q + ill, where i = we obtain: 

^ + (^ - i6C)X = {2p - l)eqJ 

Solving eq. (16) and using eq. (17) we get: 

-| 1/2 


TTi = {2p- l)e 


1 — 2e cos Tr + e 
r 2 + t 2 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


7r„ = {2p - l)^e, 


«r| + r2 


+ 


[2TrT^ cos Tr + (t^ — T^) sin — 2TrT£ 

+ t2 ^ 


^9,) 


where = S(t, Tc = Qt ^nd rj = are the rota¬ 
tion, conversion and “correlation” depths, respectively. 
As in the synchrotron thick case, circular polarization can 
exceed linear polarization. For example, when is negli¬ 
gible, the CP/LP ratio exceeds unity if: 


{2p - 1)t, > 2 


sm 



( 20 ) 


As equations (6)-(9) for the synchrotron thick case and 
equations (18) and (19) for the synchrotron thin case 
clearly demonstrate, correlations induced by the gradient 
in Faraday rotation across turbulent regions have a depo¬ 
larizing effect on linear and circular polarizations. 


4.1.3. Effect of statistical fluctuations on the mean 
Stokes parameters 

Having obtained averaged quantities, we now proceed to 
calculate the effect of statistical fluctuations on the mean 
Stokes parameters. We introduce the following notation 
for the mean and fluctuating parts of the Stokes param¬ 
eters: S = S + S, where {Si) = 0; and for the angular 
distribution of the projected magnetic field: sin2()) = a, 
cos2(/) = (2p—1)-|-&, where (a) = (6) = {ab) = 0. We write 
the fluctuating rotativity as Q = where {Cf} = 0 but 

C ~ C’(C)- We constrain ourselves to the case dominated 
by Faraday rotation, i.e., we have |CQ| ^ kCg^l) etc. In 
such a case, the lowest order terms do not depend on the 
fluctuations in the orientation of the projected magnetic 
field (j). Thus, retaining only leading terms, we obtain the 
fluctuating part of the transfer equations for Q and U: 


f + <> + S) = 0 

(21) 

dU 

(22) 

~ CiQ + Q) = 0 . 


Note that when the fluctuations in Q and U are not 
dominated by rotativity terms, they are determined by 
variations in the orientation of the magnetic field. From 
equations (21) and (22) we can obtain corrections to the 
mean Stokes parameters due to the gradient in Faraday 
rotation across each cell: 


(QQ) = - + U)'^ = -^UiCAr) (23) 

&) = (24) 

where At = AT(sin 0)“+^/2^ q jg angle between the 
line of sight and the direction of the magnetic field, At = 
Tq/N, To is the maximum synchrotron depth and N is the 
number of turbulent zones along the line of sight within 
this depth. Although these terms play an important role 
in the expressions for the mean Stokes parameters (equa¬ 
tions [6]-[10] and [18]-[19]), they introduce only higher or¬ 
der corrections ~ (!1 (C^At/ 2) to the fluctuation equations 
and therefore should not be taken into account in equa¬ 
tions (21) and (22) for the treatment to be self-consistent 
(recall that we assume that ((At < 1). Note also that both 
{QQ) and {QU) are proportional to ToA“^ and therefore 
asymptotically tend to zero as the number of field rever¬ 
sals along the line of sight increases. 

We note that the presence of inhomogeneities in the 
magnetic field can, under certain circumstances, intro¬ 
duce additional complications into the radiative transfer 
due to the tracking and coupling of plasma modes. As 
the propagating wave goes from the QL to the QT regime 
(i.e. B almost perpendicular to the line of sight) and then 
back to the QL regime, it can adiabatically adjust itself 
to the shifting nature of the eigenmodes provided that 
QvAt > 1, where Co,; = Ca^b^lnT, 111 (Bjornsson 1990; 
Thompson et al. 1994; Fuki et al. 1998), i.e., provided 
that the gradient in the magnetic field is sufficiently small. 
In the opposite limit, the propagating wave does not ‘no¬ 
tice’ any inhomogeneities. In such a case, a circularly po¬ 
larized wave will preserve its helicity as it crosses the QT 
region. Therefore, our assumption that Stokes I, Q, U 
and, in particular, V, are continuous at sharp boundaries 
between the turbulent cells, is consistent with our initial 
assumption that Co^At < I. A spatially varying magnetic 
field can also introduce mode coupling when the coupling 
constant C ~ Cog/iClv^T) ~ 1, where Qq = 2(^1 Ini'y/j^) 
(Jones and O’Dell 1977b) even in the QL region (i.e., B 
not in a direction almost perpendicular to the line of sight). 
When £ ^ 1 and propagation occurs in the QL regime, 
coupling effects are unimportant. In the opposite limit 
when C 1, and when there is no uniform component of 
the magnetic field, radiation propagates as in a vacuum. 
Although we present no formal proof, we argue that in the 
latter limit radiation will be unaffected by the fluctuating 
component of the magnetic field and will be sensitive only 
to the mean bias component of the total magnetic field. 
In a real situation, the precise value of C will depend on 
the details of magnetohydrodynamical turbulence. How¬ 
ever, since our calculations were performed assuming a 
piece-wise homogeneous medium, technically our results 
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are exact. This is consistent with claims made by Jones 
(1988), who also considered the piece-wise homogeneous 
case and performed calculations using “standard” transfer 
equations and the wave equation. He obtained practically 
identical results from the two methods even though some 
of his results formally violated the non-coupling criterion. 
In a real situation we would expect some gradient of mag¬ 
netic field across each turbulent cell but, again, in general 
the results will depend on the detailed treatment of the 
MHD turbulence which is beyond the scope of this work. 


5. RESULTS AND IMPLICATIONS 

We now consider a range of specific examples aimed 
at demonstrating predictions of our model and the con¬ 
sistency of our formulae with the results of Monte Carlo 
simulations. As mentioned in Section 3.1.1, we focus on 
cases where Faraday rotation per unit synchrotron optical 
depth Q is large. Specifically, we assume that the typ¬ 
ical mean Lorentz factor of radiating electrons 7 ~ 10^ 
and that the electron energy distribution function has a 
power-law form n(j) oc 7 “( 2 a+i)^ where a is the spec¬ 
tral index of optically thin synchrotron emission. We use 
a = 0.5 and assume that the electron distribution is cut-off 
below 7 i ~ a few. For example, for the maximum bright¬ 
ness temperature Tt, ~ lO^^A" (Readhead 1994) we have 
7 ~ ikTb/meC^ ~ 50, which corresponds to mean rota¬ 
tion and conversion per unit synchrotron optical depth of 
order ~ <5(0 ~ 3 x lO^Jlnyi/yf and Q ~ -ln( 7 / 7 *)) re¬ 
spectively, for u ~ 7^eiJ/27rmeC In order to facilitate 
comparison of analytical and numerical results we use ap¬ 
propriate synchrotron depth-weighted solid angle averages 
of coefficients of rotativity Q and convertibility Q which 
include the effect of inclination of the uniform component 
of magnetic field to the line of sight: 

(Q„(i?sin0)“+3/2)^ 
((i?sin0)“+3/2)a • ^ ^ 

Retaining only the leading terms, we get 

c; = - 2 C In + ^TT 0 - In 2 ) C ( 26 ) 

C: = |C7"^^cos0„, (27) 

2 if 

where is the angle between the line of sight and the 
direction of the uniform magnetic field component. Anal¬ 
ogously, for the correlation terms (equations [23] and [24]), 
we use 

(CQ) ^ (OT ^ ^ 2 i^ 

u Q 4 ^ 73 

5.1. Mean linear and circular polarizations 

Fig. 2 presents mean linear (upper panel) and circu¬ 
lar polarizations for the case of radiation transfer through 
a high synchrotron depth. The uppermost lines on each 
panel show analytic results for “saturated” polarizations 
corresponding to ^ < 1 for p = \ (i.e., no dispersion of 
the projected magnetic field on the sky), and filled squares 
connected by dashed lines denote results of numerical sim¬ 
ulations. Small systematic differences on Figures 2 and 3 



between the analytical curves and the numerical ones at 
higher values of 5 are due to the fact that the analytical 
formulae do not include higher order J-terms. The gen¬ 
eral behavior of these curves can be understood in simple 
terms. Linear polarization gradually decreases with S (and 
mean rotativity Q) as a result of the increasing strength 
of Faraday depolarization. Note that, contrary to common 
expectation, the mean linear polarization does not vanish 
as the result of a highly inhomogeneous magnetic field even 
though the magnitude of rotation per unit synchrotron op¬ 
tical depth is very large. As expected, circular polarization 
initially becomes stronger with an increasing component 
of the uniform magnetic field parallel to the line-of-sight. 
However, as Faraday depolarization gradually eliminates 
linear polarization (both Stokes Q and U are affected by 
this process), there is a reduced amount of U available for 
conversion into circular polarization (Stokes V). Thus cir¬ 
cular polarization has an extremum. 

Other curves on Fig. 2 represent the case of “unsatu¬ 
rated” polarization when ^ ^ 1. In such situations, the 
levels of both linear and circular polarizations systemat¬ 
ically decline as N decreases. Note that the number of 
turbulent zones along the line of sight cannot be arbitrar¬ 
ily low, as this would reduce the mean polarization levels 
to very low values due to the increasing influence of the 
correlations between rotativity and Stokes Q and U pa¬ 
rameters. Apart from affecting the mean polarization lev¬ 
els, increasing the number of field reversals along the line 
of sight also reduces the polarization fluctuations. Due to 
the finite size of the telescope beam, the fluctuations are 
further suppressed, however the beam averaging process 
does not affect the mean level of circular polarization. It 
also does not influence the mean level of linear polariza¬ 
tion, provided that LP fluctuations are not very large prior 
to the beam averaging. Assuming that the number of tur¬ 
bulent zones in each of the perpendicular directions across 
the telescope beam is very roughly comparable to N, we 
conclude that fluctuations in the mean levels of CP and LP 
due to the stochastic nature of the plasma should be rela¬ 
tively small. Indeed, even though our numerical results for 
the mean circular and linear polarizations were obtained 
for a number of “pencil beams” much smaller than N'^, 
none of our results exhibit strong fluctuations. Thus, by 
demanding that there be enough turbulent cells along the 
line of sight to avoid correlation depolarization, we also 
guarantee that the polarization variability should be dom¬ 
inated by changes in the mean plasma parameters (e.g., 
the magnitude of the uniform component of the magnetic 
field, dispersion of the projected random magnetic field 
on the sky, the synchrotron depth, etc.), rather than by 
statistical fluctuations in, e.g., the orientation of the local 
turbulent magnetic field. This also assures that the sign 
of circular polarization should be a persistent feature — 
variations of the mean plasma parameters are not likely 
to change the sense of circular polarization, provided that 
the synchrotron optical depth does not change dramati¬ 
cally from low to very high values (see below) and the 
orientation of the mean field remains the same. 

Fig. 3 illustrates the effect of relaxing the assumption of 


^ We further comment on the choice of parameters and consider more general situations while discussing specific observational cases in Sections 
6.1 and 6.2 
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that the the projected magnetic field has zero dispersion 
on the sky. As the projected orientation of the magnetic 
field becomes increasingly chaotic {p —>■ 0.5), the magni¬ 
tudes of linear and circular polarizations gradually decline. 

As both Figures 2 and 3 indicate, circular polarization 
can exceed linear polarization at higher values of the bias 
parameter 6. The excess of CP over LP requires signifi¬ 
cant Faraday depolarization of linear polarization. Never¬ 
theless, even a small amount of U can then be effectively 
converted to circular polarization, leading to CP/LP ra¬ 
tios in excess of unity. It is even conceivable to have a 
situation in which linear polarization falls below the detec¬ 
tion threshold whereas circular polarization is still easily 
observable. Note that significant Faraday depolarization 
does not require a large bias S (i.e., exceeding unity) be¬ 
cause of the large values of rotativity. This is consistent 
with observations which do not reveal any dominant large 
scale unidirectional fields. Such fields would lead to very 
strong linear polarization and are also unexpected on the¬ 
oretical grounds. 

In real sources the effective synchrotron optical depth 
will depend not only on the emission properties of the 
plasma along the line of sight but also, among other fac¬ 
tors, on the solid angle subtended by the emitting region. 
In a realistic situation, the telescope will integrate over a 
finite-size beam with different synchrotron depths along 
different lines of sight. Thus, the emission from the syn¬ 
chrotron self-absorbed core of a jet will be weighted by 
a smaller solid angle than the emission from the more 
extended regions which have lower synchrotron depths 
and, therefore, lower surface brightness. Detailed calcu¬ 
lations of the polarization properties of specific jet models 
are beyond the scope of this paper and will be presented 
in a forthcoming publication (Ruszkowski and Begelman, 
in preparation). Fig. 4 illustrates the effect of varying 
the total synchrotron depth Tq for ‘saturated’ polariza¬ 
tion < I). For large Tq the results are similar to the 
ones presented on Fig. 2 and 3. In the opposite extreme, 
i.e., when Tq tends to zero, we have qualitatively similar 
behavior with the main difference being that the helicity 
of circular polarization is reversed. As in the high syn¬ 
chrotron depth case, the general trends in the polarization 
behavior can be understood in simple terms. The grad¬ 
ual decline of linear polarization with 6 is just a result of 
Faraday depolarization. As radiation propagates through 
the plasma, Stokes U is generated from Stokes Q due to 
Faraday rotation. Therefore, circular polarization, which 
is produced from Stokes U by Faraday conversion, initially 
increases. However, for large S, Faraday depolarization re¬ 
duces Stokes Q and U and thus leads to the suppression 
of circular polarization. The oscillations in linear and cir¬ 
cular polarizations are due to cyclic rotation of Stokes Q 
into U followed by conversion to V. Note that, for a nar¬ 
row range of intermediate values of synchrotron depth, the 
behavior of CP resembles that of low Tq when S (and the 
mean rotativity) is small, and that of high Tq when S is 
larger. 

The effects of “correlation depolarization” for syn¬ 
chrotron depths around unity are showed in Fig. 5. This 
figure also demonstrates our analytic results for small syn¬ 
chrotron depths (dashed lines). Clearly, our analytic so¬ 
lutions for both very large and small synchrotron depths 


adequately describe the numerical results which smoothly 
join the two regimes. 

6. APPLICATION TO SPECIFIC SOURCES 
6.1. The case of quasar 3C279 

Wardle et al. (1998) reported the discovery of circu¬ 
lar and linear polarization in 3C 279 and attributed CP 
to internal Faraday conversion. Typical fractional linear 
and circular polarizations in 3C 279 are of order ~ 10% 
and < 1%, respectively. Wardle et al. (1998) concluded 
that if the jet is composed of normal plasma, then the 
low-energy cut-off of the energy distribution of relativistic 
electrons must be as high as ~ 100 in order to avoid 
Faraday depolarization and overproduction of the jet ki¬ 
netic power. They considered synchrotron thin models 
with a bias in magnetic field and a small number of field 
reversals along the line of sight but neglected the disper¬ 
sion of the projected magnetic field on the sky. They were 
unable to fit their polarization models to the observational 
data for > 20 and thus claimed that the jet must be pair- 
dominated. However, the above observational constraints 
on CP and LP and the jet energetics can be satisfied for 
a variety of microscopic plasma parameters. This is due 
to the fact that different “microscopic” parameters, such 
as 7 i, the ratio of the cold to relativistic electron num¬ 
ber densities, or the positron fraction, can lead to similar 
“macroscopic” parameters such as convertibility and rota¬ 
tivity. In order to illustrate this, we consider two radically 
different examples and show that both cases can lead to 
the same CP and LP. 

6.1.1. Electron-proton jet 

In this example, plasma is composed exclusively of a 
mixture of protons and electrons with both relativistic 
and cold populations being present. For instance, for 
a low-energy cut-off ~ 30 and an electron number 
density-weighted mean Lorentz factor y ~ 50, we get 
iCv) ~ 80{nc/nr)S and (Q) ~ —0.5 (for Uc/ur <C 10^; 
see Appendix for details). The above value of rest-frame 
y is consistent with self-absorbed sources having bright¬ 
ness temperatures in a narrow range close to ~ lO^^iL 
(Readhead 1994). The actual energy distribution of the 
radiating particles will not be characterized by sharp en¬ 
ergy cut-offs but will rather be a smooth function (Shi- 
mada and Hoshino 2000; Dieckmann et al. 2000). How¬ 
ever, the detailed treatment of these subtleties is beyond 
the scope of this paper and we believe that the overall 
complexity of the problem justifies the use of our approx¬ 
imate treatment. For the above choice of parameters, the 
main contribution to rotativity comes from cold electrons 
as long as Ucfur ^ 5 x 10“^ (cf. eq. [A9]). The required 
levels of LP and CP can be obtained, for example, when 
p = 1 and the combination {nc/nr)S ~ 2.5 x 10“^, as 
this assures sufficiently low mean rotativity that it does 
not lead to Faraday depolarization. Bear in mind that 
the admixture of cold electrons does not have to be large 
to explain the data. For example, we get the right levels 
of LP and CP for Uc/ur ~ 5 x 10“^ and (5 ~ 5 x 10“^. 
Interestingly, a jet with such a plasma composition could 
carry roughly as small a kinetic power as the pure electron- 
positron jet with the same emissivity, since the ratio of 
kinetic powers of an e-p jet to a pure relativistic 
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jet is ~ 18.4((7)e+e-/50) ^( 7 i,e+e-/ 7 i,pe), where we have 
assumed that protons are cold and have used a = 0.5 
and where e“'"e“ and pe refer to pair plasma and normal 
plasma, respectively. This means that, in principle, one 
can also explain the mean levels of circular and linear po¬ 
larizations in 3C 279 for much lower values of 7 ^ than the 
one used above and therefore much higher rotativity, pro¬ 
vided that the bias parameter 5 is also much lower. For 
example, the right levels of CP and LP can be obtained 
for 7 j = 3, 7 = 50, = 1875, (5 ~ 3 x 10“^ and t = 1 (see 

middle panels in Fig. 5). We note that observational con¬ 
straints on jet composition based on energetics should be 
taken with caution in any case, as kinetic luminosities de¬ 
rived from these methods scale with high powers of poorly 
determined observational quantities. 

6.1.2. Electron-positron jet 

The alternative possibility is that the jet is dominated 
by relativistic pair plasma. For example, for 7 ^ = 2 and 
7 = 50 we get (C*) ~ -3.1 and (Q) ~ 1.4 x 10‘^S{np/ne), 
where Up and Ug are the number densities of protons and 
electrons, respectively. Agreement with the observed frac¬ 
tional linear and circular polarizations can be obtained for 
p = 1 and {np/ne)6 ~ 2.1 x 10“^. Depending on the 
actual values of S and the ratio of protons to electrons, 
the jet may be pair-dominated in the sense that rig Up 
while being dominated dynamically by protons; or it can 
be dominated by pairs both numberwise and dynamically. 
Recent theoretical work of Sikora and Madejski (2000) 
suggests that jets may be pair-dominated numberwise but 
still dynamically dominated by protons. 

6.2. The case of Sgr A* 

Among the most intriguing sets of polarization observa¬ 
tions are those of the Galactic center. Observations of stel¬ 
lar proper motions in the vicinity of the nonthermal source 
in the center of the Galaxy (Sgr A*) reveal the presence 
of a ~ 2.6 X 10®Mq compact object — the most convinc¬ 
ing candidate for a supermassive black hole (Eckart and 
Genzel 1996, 1997; Ghez et al. 1998). Recently Bower et 
al. (1999) reported the detection of circular polarization 
from Sgr A* with the VLA, which was confirmed by Sault 
and Macquart (1999) using ATGA. The typical level of 
GP in their observations was ~ 0.3%, greater than the 
level of linear polarization. This result may seem sur¬ 
prising in light of the strong limits on the ratio of GP to 
LP in AGN where GP/LP is usually much less than unity 
(Weiler and de Pater 1983). However, as explained above, 
an excess of CP over LP can be explained easily in the 
framework of our model. Archival VLA data indicate that 
the mean CP was stable over ten years (Bower 2000). 
This is also not surprising as our model naturally predicts 
a persistent CP sign provided that the number of field 
reversals, either along one line-of-sight or across the beam 
area, is sufficiently large and the source does not undergo 
dramatic changes from synchrotron thin to synchrotron 
thick regimes. However, short-term circular polarization 
variability (~ a few days) is also present. Interestingly, 
recent VLA observations from 1.4 to 15 GHz separated 
by a week revealed a GP increase at frequencies greater 
than 5 GHz, which coincided with an increase in the total 
intensity (Bower et al. 2000). The CP spectrum was 


characterized by a flat to slightly positive spectral index 
(tTc oc > 0). This result can also be accounted for 

in our model. For example, in the framework of a self- 
absorbed, self-similar jet model (Blandford and Konigl 
1979) but with a small bias S we have {Q) oc S oc B oc v. 
Assuming that ^ > eqfq and q = {2p — l)^(Cg)^ 1, and 

then demanding ItTcI/tt/ > 1 and (3 0, we obtain the 

following criterion which leads to the observed behavior 
(cf. eq. [7]-[10]): 

(C + < (<50" < 2C + <? + 0? + <?)■ (29) 

For instance, for ^ = 20, (Q) ~ —10, S( ~ 40 and p = 1 
we have /3 > 0, {tTcI/tti ~ 3 and \tTc\ ~ 0.6%. Such values 
of the mean convertibility and rotativity could be obtained 
in the case of predominantly cold e~-p plasma with a small 
admixture of relativistic electrons. Note, however, that a 
certain amount of relativistic particles is always essential 
as synchrotron processes dominate emission and absorp¬ 
tion in typical nonthermal sources. Specifically, the min¬ 
imum and maximum size of Sgr A* constrain the bright¬ 
ness temperature to be 10^® ^ Tj, < 5 x 10^^AT (Melia 
and Falcke 2001), which is within the range of typical 
AGN radio cores. Taking ~ 10^^ as the representative 
rest frame value (Readhead 1994), we get 7 ~ 50. For 
minimum cut-off Lorentz factor 7 j = 2, this gives rela¬ 
tivistic contributions to convertibility and rotativity equal 
to (Cg^^^) ~ —3 and ~ 2755, respectively. Pro¬ 

vided that dielectric suppression (Razin effect) is unim¬ 
portant, the overall transfer coefficients are merely the 
sums of cold (c) and relativistic (r) contributions. Thus, 

for ~ 4.3(nc/nr) and ~ 0.08(nc/nr) 

(see Appendix), we obtain (nc/rir) ~ 30 and 5 ~ 10“®, 
which corresponds to the required values of the mean con¬ 
vertibility and rotativity. 

It has recently been suggested that observations of lin¬ 
ear polarization can be used to constrain the accretion 
rate in Sgr A* and other low-luminosity AGN (Agol 2000; 
Quataert and Gruzinov 2000). These authors base their 
argument on the assumption that the Faraday rotation 
measure has to be sufficiently small in order not to sup¬ 
press strong linear polarization (Aitken et al. 2000). This 
assumption places limits on density and magnetic field 
strength and leads to very low accretion rates ~ 10 “® to 
IO“®M 0 yr“^. As noted by Agol (2000) and Quataert and 
Gruzinov (2000), this is inconsistent with an advection- 
dominated model for Sgr A*, which assumes that the ac¬ 
cretion rate is of order the canonical Bondi rate ~ 10“^ to 
\0~^Mq yr“^. We point out that strong rotation measure 
does not in principle limit densities and magnetic fields 
provided that the field has a very small bias superimposed 
on it, which is required to define the sign of circular po¬ 
larization, and that the dispersion of the projected mag¬ 
netic field on the sky is not too large. This implies that 
‘high’ accretion rates comparable to the Bondi rate can¬ 
not be excluded on these grounds. More work on this 
issue is required to fully exploit the information contained 
in polarization observations in order to constrain physical 
conditions in Sgr A*. 

6.3. Radio galaxy M8T 
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Brunthaler et al. (2001) detected circular polarization 
in the compact radio jet of the nearby spiral galaxy M81. 
Their estimated values of CP were 0.27 ± 0.06 ± 0.07% at 
4.8 GHz and 0.54 ± 0.06 ± 0.07% at 8.4 GHz, where errors 
are separated into statistical and systematic terms. This 
suggests that the CP spectrum is flat or possibly inverted. 
They also detected no linear polarization at a level of 0.1%, 
indicating that the source has a high circular-to-linear po¬ 
larization ratio. The spectral index indicates that this 
source is synchrotron thick and we can apply the same ap¬ 
proach as for Sgr A*. For example, for ^ ~ 50 and ~ 80 
we get a flat to slightly positive CP spectral index, CP/LP 
ratio ~ 4.7 and tTc ~ 0.33%. 

6.4. IDV hlazar PKS 1519-273 

Macquart et al. (2000) reported strongly variable cir¬ 
cular and linear polarization from the intraday variable 
blazar PKS 1519-273. The source was characterized by 
an inverted intensity spectrum. The spectrum of circular 
polarization was roughly consistent with being flat or in¬ 
verted (see their Fig. 3c) and typical values of linear and 
circular polarizations were a few percent and ~ 1%, respec¬ 
tively. These data can be explained, for example, in the 
limit of high synchrotron depth. In order to have a roughly 
flat/inverted CP spectrum, we need the bias to correspond 
approximately to 6 near or below the extremum in circular 
polarization according to curves in Fig. 3. Levels of linear 
and circular polarization, in qualitative agreement with 
the observed ones, can be obtained for various combina¬ 
tions of ^ and the dispersion of the projected random mag¬ 
netic field on the sky (2p — 1) (cf. equations [6]-[8] and see 
Figures 2 and 3). Alternatively, the data can be explained 
for synchrotron depth r ~ 1. In this case, the flatness 
of the synchrotron spectrum could be explained by means 
of a Blandford-Konigl model in which the photospheric 
surface of constant, model-dependent, synchrotron depth 
moves outward along the jet with decreasing frequency. 
The required levels of CP and LP can then be obtained 
easily for the parameters considered in Fig. 5 (see middle 
and right panel). 

6.5. X-ray binary SS 433 

Fender et al. (2000) detected circular polarization from 
the radio jet in the famous X-ray binary SS 433. The 
flux density spectra of circular polarization and Stokes / 
were roughly of the form V (x I (x . However, mul¬ 
tiple components in the source and a lack of high spatial 
resolution prevented determination of the origin of circu¬ 
lar polarization and the spectrum of fractional polariza¬ 
tion. They argued that the CP emission is likely to be 
produced in the innermost regions of the binary, and that 
the fractional CP of this region can be as high as 10%. In 
this case the ratio of the fractional circular polarization to 
the overall linear polarization (< 1%) would exceed unity. 
If circular polarization originates from regions character¬ 
ized by synchrotron depth t ^ 1, then the underlying 
synchrotron spectrum from this region would be approx¬ 
imately flat/inverted (Paragi et al. 1999). Then, the 
fractional CP spectrum would scale roughly oc which 
corresponds to a bias 6 (and mean rotativity) greater than 
that associated with the CP extremum. The highest cir¬ 
cular polarization would then be about 3%. Note that in 


such a situation we would also expect an excess of CP over 
LP even in the innermost regions of SS 433. If CP origi¬ 
nates from the inner parts of the binary but from regions of 
r < 1, then the flatness of the synchrotron spectrum could 
be explained by means of a Blandford-Konigl model. In 
such a scenario, the maximum fractional CP from the cen¬ 
tral regions could be as high as ~ 20% (see Fig 4) and the 
CP spectrum could have a negative slope [3 (tTc oc cf. 
eq. [19] for rj <C 1). 

7. CONCLUSIONS 

We have considered the transfer of polarized syn¬ 
chrotron radiation in stochastic sources by means of an 
analytic approach and a set of numerical simulations, and 
have argued that Faraday conversion is the primary mech¬ 
anism responsible for the circular polarization properties of 
compact radio sources. A crucial ingredient of our model 
is a small bias in the highly turbulent magnetic field which 
accounts for the persistence of the sign of circular polar¬ 
ization. This bias is direct evidence for the net magnetic 
flux carried by magnetically accelerated jets (e.g., Bland- 
ford and Payne (1982); Li et al. (1992)). 

Extremely large rates of Faraday rotation, i.e., Fara¬ 
day rotation per unit synchrotron absorption depth, do 
not necessarily lead to depolarization provided that the 
mean rate of Faraday rotation across the source is rela¬ 
tively small, or in other words, that the turbulent mag¬ 
netic field possesses a very small directional bias. Indeed, 
a large Faraday rotativity is required in order to explain 
the high ratio of circular to linear polarization observed in 
some sources. Gonstraints on jet composition or accretion 
rate, based on the requirement that the source does not 
become Faraday depolarized, may be circumvented under 
these conditions. 

Gradients in Faraday rotation across turbulent cells can 
lead to correlations between rotativity and Stokes Q and 
U parameters, which can result in “correlation depolariza¬ 
tion”. Observed polarization levels require that the field 
have many reversals along the line of sight to avoid this 
effect. Statistical fluctuations of circular and linear polar¬ 
izations are then likely to be dominated by changes in the 
mean parameters describing the plasma rather than by the 
stochastic behavior of the turbulent medium. Variations 
in the mean parameters are unlikely to change the helicity 
of circular polarization unless a source undergoes a sharp 
transition from very low to very high synchrotron depth. 

We have shown that our model is potentially applica¬ 
ble to a wide range of compact synchrotron sources. In 
particular, it naturally predicts an excess of circular over 
linear polarization when a source is strongly depolarized 
by the mean Faraday rotation and when a small amount 
of linear polarizarion is efficiently converted into circular 
polarization. This can explain the polarization properties 
of the Galactic Center and M81*. 

We thank Roger Blandford, Avery Broderick and Marek 
Sikora for insightful discussions. This work was supported 
in part by NSF grant AST-9876887. 



Circular polarization from synchrotron sources 9 

APPENDIX 

TRANSFER OF POLARIZED RADIATION 

The transfer equation of polarized radiation reads (Zheleznyakov et al. 1974; Jones 1988) 

/ (^ + 1) C<?cos2(/) -(gsm2(j) \ / I \ 

C, cos 20 (^ + 1) C Q sin 20 Q 

-C, sin20 -C i^ + l) Q cos 20 U 
V C. -q sin 20 -q cos 20 + 

where I, Q, U and V are the usual Stokes parameters, r oc is the synchrotron optical depth, J is the 

source function, 0 is the azimuthal projection angle of the magnetic field on the sky and the coefficients of emissivity (e^, 
Cy), absorptivity (q, q), convertibility (q) and rotativity (0*) are given by: 






(A2) 

(A3) 

(A4) 

(A5) 

(A6) 


(A7) 


where C^q /C%'^ 
by: 


are the ratios of convertibilities {q) and rotativities (v) of the cold (c) and relativistic (r) plasmas given 


^ _ J_ g- 1/2 j_ n- +n+ 

^*{r) 2a 1 - (p/Pi)“-l/2 ry. 

^ J_ Q: + 1 jI_ n- - n+ 
q(’') 2a a + 3/2 Iny^ Ur — nt 


(AS) 

(A9) 


The transfer coefficients have been generalized to include contributions from both electron and positron plasmas but 
they assume an isotropic pitch-angle distribution of the radiating particles. We also neglected circular emissivity and 
absorptivity (y (Jones and O’Dell 1977a). In the above equations, n+q are number densities of cold/relativistic electrons 
(—) or positrons (-I-), i^b = eB sin 0/27rmeC is the Larmor frequency, 7^ is the low-energy cut-off Lorentz factor of relativistic 
particles, Pi = y/pe is the frequency corresponding to radiating particles of energy ~ jirUeC^, a is the synchrotron-thin 
spectral index which also defines the slope of the relativistic particle energy distribution 71(7) oc 7“^““^, 9 is the angle 
between the line of sight and the direction of the magnetic field and e® , q, 0*® and 0*” are the proportionality coefficients 
which are tabulated in Jones and O’Dell (1977a). 

We integrated the radiative transfer equations using the Cash-Karp embedded Runge-Kutta fifth-order method with 
an adaptive stepsize control (Press et al. 1992). Computations were performed by integrating transfer equations in a 
piece-wise homogeneous medium. The physical conditions in each cell, i.e., orientation of the magnetic field vector, were 
chosen using the random number generation method of L’Ecuyer in the implementation provided by Press et al. (1992). 
This method is particularly suited for our purposes as it generates long-period sequences of random numbers and thus 
prevents any spurious correlations. We assumed a constant strength of the variable B-field component but allowed its solid 
angle distribution to be uniform within 9 G (0°, 180°) and 0 G (— 0o,0o), where (cos 20)(_0^ = 2p — 1 is the average 

of cos 20 in the interval (— 0o,0o). We also superimposed a weak constant and unidirectional component |B„| = (5|B| on 
the variable magnetic field. 


® Our formulae are free from typographic errors found in Jones (1988). 
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Fig. A2.— Linear (upper panel) and circular polarizations for high synchrotron depth (tq ~ 27) and varying number of turbulent zones 
along the line of sight {N = 10®, 3 X 10® and 10®; from top to bottom, respectively). Other parameters used in this simulation were: a = 0.5, 
7 i = 3, 7 = 50, p = 1 and = 45°. Solid lines denote analytical results and squares connected by dashed lines denote results of Monte Carlo 
simulations. 
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Fig. A3.— Linear (upper panel) and circular polarizations for high synchrotron depth, large number of field reversals along the line of sight 
(‘saturated’ polarization) and varying dispersion of the projected turbulent magnetic field on the plane of the sky (</> is distributed uniformly 
in [—A(/>, A((>]; A(/> = 0“’ (p = l)i = 30° {p = 0.91), A<f> = 45° (p = 0.82), A(f> = 60° (p = 0.71); from top to bottom, respectively. Other 

parameters used in this simulation were: a = 0.5, 7 i = 3, 7 = 50 and = 45°. Solid lines denote analytical results and squares connected 
by dashed lines denote results of Monte Carlo simulations. 








Circular polarization from synchrotron sources 


15 



T= 1 .0 


t=2.0 



0.01 


0.1 0 


0.01 


0.1 0 


0.01 


0.1 0 


Fig. A5.— Linear (upper panels) and circular polarizations (lower panels) for moderate values of synchrotron depth (top labels) and 
different numbers of turbulent cells N along the line of sight. The curves correspond to N = 187, N = 562 and N = 1875 (left panel; from 
bottom to top, respectively), N = 625, N = 1875 and N = 6250 (middle panel; from bottom to top, respectively) and N = 1250, N = 3750 
and N = 12500 (right panel; from bottom to top, respectively). Other parameters used in this simulation were: a = 0.5, = 3, 7 = 50, 

p = 1 and Ou = 45°. Dashed lines on the left panel are the analytical results for small synchrotron depth. 







